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Abstract. - The question as to whether or not quantum mechanics is applicable to the macroscopic 
scale has motivated efforts to generate superposition states of macroscopic numbers of particles and 
to determine their effective size. Superpositions of circulating current states in flux qubits constitute 
candidate states that have been argued to be at least mesoscopic. We present a microscopic analysis 
that reveals the number of electrons participating in these superpositions to be surprisingly but not 
^ ' trivially small, even though differences in macroscopic observables are large, 
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Introduction. — Quantum mechanics makes predictions that often seem at odds with 
the intuition gained through human experience at the macroscopic level. Schrodinger's 
cat paradox pQ provides an extreme example of this, highlighting the consequences of ex- 
trapolating basic quantum concepts such as superposition from microscopic to macroscopic 
systems. The question of whether quantum behavior is restricted for large numbers of par- 
ticles by some unknown non-quantum mechanism or contains some limitation that we do 
not yet understand is a fundamental unresolved question [2]. Realizing true macroscopic 
quantum superpositions would constitute one step towards an answer, providing evidence 
against macroscopic realism [3J. 

There are many different ways in which one could define the effective "size" of a macro- 
scopic superposition state. One might first compare the absolute value of the difference in 
some suitable physical observable between the two branches to some characteristic atomic 
scale for that observable. However, as discussed in Ref. [2], this measure is not an unam- 
biguous signature of a large quantum superposition of microscopic degrees of freedom. In 
order to understand how many microscopic degrees of freedom actively participate in the 
superposition, it is essential to further make a microscopic analysis of the differences be- 
tween the two branches of a superposition state. In this work we provide for the first time 
such a microscopic analysis for the superposition states made with flux qubits that show 
macroscopically distinct currents in different branches and which have consequently been 
suggested to involve superposition of macroscopic numbers of electrons [I] . 
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Superconducting flux qubits composed of micron-scale superconducting loops containing 
Joscphson junctions and involving currents attributable to f0 6 -10 10 electrons 05H7] provide 
good candidates for forming macroscopic quantum superpositions. Recent experiments have 
demonstrated that flux qubits can be brought into a superposition of macroscopically distinct 
states of circulating currents, of the form |0) + |0) [HE], and that coherent oscillations 
between these states can be induced [8]. Here, the states |0) and |0) correspond to states 
of fluctuating persistent current with nonzero mean values in clockwise and anticlockwise 
directions around the loop, respectively [9]. However, while the two components of this 
superposition state have macroscopically distinct values of magnetic flux and operation of 
flux qubits may be described in terms of a few macroscopic variables [5H51[Tfj]. it is not 
clear how many microscopic degrees of freedom participate in the superposition [2ll4lll 1 j. 
Answering this fundamental question requires a fully microscopic analysis for the number of 
electrons participating in the superposition, a challenging task which has hitherto not been 
undertaken for superconductors. We present here a first principles approach that allows 
us to place a bound on this size and evaluate this for experimental flux qubit systems. 
We find that for all flux qubits fabricated to date, this results in a superposition state 
whose size is considerably smaller than the number of electrons carrying the supercurrent. 
Our microscopic analysis reveals that the indistinguishability and fermion statistics of the 
constituent Cooper pair electrons play a critical role in this reduction. 

Our calculations are made at T = and are based on the framework of BCS theory [12] . 
in which pairing is described by a complex- valued order parameter A(r) = |A|e^. In order 
to describe a flux qubit, we use the path integral description |13j of a BCS superconductor 
in a ring geometry, interrupted by small capacitance Josephson junctions. In this geometry, 
Coulomb repulsion at the junction drives quantum fluctuations of 4>{r). Cat states in flux 
qubits are described as superpositions of many-body states with clockwise and counterclock- 
wise net circulating currents, which we denote by A= + |t>))- In each of the two branches 
of the superposition, the phase is given by a narrow distribution with mean value </>L/ij(r), 
respectively. For qubits with thickness much smaller than the penetration depth A, we can 
assume that the current flow is uniform in the lateral dimension (the thickness, or vertical 
dimension of a planar geometry loop with horizontal orientation). The order parameter and 
the vector potential A together define the superfluid velocity vs = — ^7 (V</> + i?A), with 
m the single electron mass and e the single electron charge, vs is equal to the mean velocity 
of superconducting electrons in the system as normally defined in terms of the center of 
mass velocity of Cooper pairs and is related to the electron current density via j = ep e vs, 
with p e an effective electron density [14] . (The Cooper pair superfluid density is equal to 

Pell) 

As long as both the externally applied and the internally generated magnetic fields 
are weak enough and the cooling is sufficiently adiabatic to avoid vortex formation, we 
can decompose the entire electron system in terms of single-electron modes labeled by the 
linear momentum p = ftk and spin er, where k is the eigenvalue of — iV. To analyze 
the current carrying states, we define creation and annihilation operators CT and c q (T , 
where q is the eigenvalue of — iV + 7^ vs. The index k labels the internal momentum of 
Cooper pairs. In BCS theory, the average occupation number of a single-electron mode is 

given by = (d^ ^.Ck^y (rik is independent of a for any realistic magnetic field strength) 

and at zero supercurrent is equal to rik = (1 — -Efc/f2fc)/2 with E^ — ^M: fip, fik = 

y/ + A 2 and /if the Fermi energy, while the amplitude of the Cooper pair mode (the 
'condensation amplitude' [T5]) = (^k -|.cl k A has a smooth peak of width ~ A/vp at 

the Fermi momentum fc^ [15] . We emphasize that the circulating current is a dissipationless 
supercurrent. The modes k and q should not be confused with quasiparticle states that 
carry dissipative processes and that are separated from the condensate by an energy gap. 
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Microscopic Analysis of Flux Qubit Superposition Size. We make here a mi- 
croscopic, all-electron analysis of the superposition size using the distinguishability measure 
of Rcf. [TB] . The latter provides an operational measure that asks what is the largest num- 
ber of elementary constituents in an TV particle system that need to be measured, within a 
given precision, to collapse the superposition to a single branch. Ref. |16j used this to show 
that the number of constituents that participate actively in the superposition, i.e., that do 
something different in the two branches, is given by N/n, where n is the smallest number 
of particles that must be measured to guarantee successfully distinguishing between the two 
branches. For flux qubit states, the electron indistinguishability must additionally be taken 
into account, which is done by extending the distinguishable particle analysis of Ref. |16j to 
measurements of electron mode occupation, with the density matrices replaced by fermionic 
multi-mode density matrices [17] . The microscopic number of constituents (modes) actively 
participating in the superposition is then determined by the number of n-electron modes 
that are required to be measured in the N mode electron system to successfully distinguish 
the two branches. Expanding the reduced mode density matrices in the plane wave basis 
q, we find that for small differences in mode occupancy between L and R branches, the 
limiting number of mode measurements can be bounded by N/AN to t where AiV to t is the 
total difference in mode occupation number between the two branches [18], leading to a 
bound AiVtot on the superposition size. This simple-looking result is underscored by a more 
profound fact, namely that when the indistinguishability of electrons is taken into account, 
ANtot is the only meaningful indicator of how many electrons are affected when passing 
from one branch of the superposition to the other. 

It can be shown that this choice of single-particle basis maximizes the occupation number 
differences between the two branches, <5n q = (0| c qcr Cq iCr O) — (O c qo .c qj(T |0) [E], thereby 
ensuring that AiV to t is a mode-independent and hence true upper bound on the cat size. 
We show below that this calculation of superposition size in terms of single electron modes 
can be directly related to the corresponding calculation in terms of Cooper pair modes. 
While the latter appears more natural, it requires more complex calculations and, as we 
discuss below, for small superfluid velocities vs, gives identical results to calculations based 
on single electron modes. 

ANtot can be calculated microscopically in the path integral approach within the BCS 
description. The average occupation number of a single-electron mode is given by ?i q = 
(c q cr Cq iCr ) (as noted above, n q is independent of a for any realistic magnetic field strength). 
The necessary expectation values can be obtained from the generating functional 

Z[$ = J ^DAD^e^ 1 *^* 1 ^^ T / ^--[et^Wt^t+fctrTlCtr.T^ ( X ) 

where VKV^) are the electron field operators, A the vector potential, ip the electromagnetic 
scalar potential and S[ip, A, <p] the action functional of Ref. [13]. Employing a Hubbard- 
Stratonovich transformation and saddle point analysis as in that work introduces the order 
parameter A(r) and reduces the generating functional to a form in which only the phase <f> 
remains as a non-classical variable, namely, 

Z[£] = y , p0 e -^^[5 C * 2 -^ cos W+2t(*-*oxt) 2 +/d 3 rt l G^ (2) 

Here the flux $ is related to the phase variable (f> through the flux quantization relation 
$ = y~ $o , the Josephson energy Ej is an effective quantity [TO] and G is the sum of the 
zcroth order bulk Green's function and perturbativc contributions from the superfluid flow 
and junction tunneling, G = (G^ lk + T rr > ) = Go + 5G Vs + SGt- Standard analysis 
then leads to the single electron and Cooper pair mode correlation functions in terms of 
matrix elements of G [20]. We analyze first the loop current contribution Go + SG VS and 
discuss the tunneling contribution SGt separately below. The loop current contribution to 
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the single mode correlation functions is obtained as 

(*|c qCT c qCT |*) = -h lim G^q-^O^-^.T), (3) 



where G aa is a matrix element of Go + <5G VS indexed by the electron spin. As long as the 
current in the superconductor is small compared to the critical current I c .buik, the supcrfluid 
velocity vg will be a small perturbativc quantity and the Green's functions can be expanded 
to first order in vs|/w C rit = vs|ufto/A, where v-p is the Fermi velocity. Carrying this out 
yields PUDS] 

1 A E q \ 1 A 2 



W^|*) = -^l-^j + -_ ftq .(v s ), (4) 

where (vg) is the mean supcrfluid velocity averaged over the quantum state of the system. 
To first order in vs, the difference in occupation number of a mode (q, a) between the two 
circulating current superposition branches, is then given by 

A 2 

<K = (01 c q , ff c q , ff |0) - (C5| c q , CT c q , CT |0) = • ^(vs) , (5) 

where 5{vs) is the difference in supcrfluid velocity between the two branches. 
Analysis of the corresponding Cooper pair correlation functions 

AA/k,-k = (0| 4t e - ki c-kic kt |c3) - (0| cj^cl^c-k+ckt |0) (6) 

in the laboratory frame is complicated by the fact that the modes q = k + mvg/fi and 
q' = k' + mvs/fo will be coupled only if k' = k, which means that whether two modes 
q and q' are Cooper pair coupled or not depends on the supcrfluid velocity. Detailed 
calculation shows that if the modes are not Cooper pair correlated, there is no first order 
contribution to the correlation function <5A/q. q ' , while if they are Cooper pair correlated, the 

corresponding occupation number A/k = ^C^CjA = ^c^ ^cL k ^c_k,j,Ck,-|-\ of a Cooper-pair 
mode is identical to [17j . Hence, looking at Cooper-pair modes rather than single- 
electron modes does not change the effective cat size and from now on we shall consider just 
the single-electron quantity <5n q , eq. ([5]). 

Superposition Size Estimates. — The change of the mode occupation Sriq reflects 
the fact that a supercurrent gives an additional momentum 2q to each Cooper pair, thereby 
changing the momentum structure of the condensate without exciting quasiparticles. To 
obtain numerical estimates, we proceed by first determining the number of electrons changing 
modes per unit volume, 5n(r), at a point r in the superconductor, by integrating the absolute 
value of <5n q over all q values. (At T = 0, there is no contribution from quasiparticle 
excitations.) We can similarly obtain an expression for <5j(r), the difference in local current 
density, as an integral over q. These two integrals can be simplified by noting that the 
contribution of modes far from the Fermi surface is strongly suppressed in eq. ((4]) and hence 
also in <5n q . As shown in Ref. |T5] , this leads to the relation Sn(r) = 3 |<5(j(r)) /4 evp, where 
we have divided by 2 in order to avoid double-counting electrons when they are removed 
from one mode and added to another. Integrating this expression over the entire volume of 
the superconductor then gives the total difference in occupation numbers as 

3L 

AJV to t = ~ A Sip, (7) 

4 evF 

where L is the total length of the main superconducting loop of the flux qubit and SI P the 
difference in "persistent current" [5] between the superposed branches. The experimental 
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Table 1: Parameters and effective cat sizes for current superposition states produced at SUNY [7J, 
Delft [B] and Berkeley [5T|. vp is the Fermi velocity, L the length of the superconducting loop, 
5I P the measured difference in persistent current between the two branches and 5/j, — ASI P is the 
difference in magnetic moment, where A is the area enclosed by the loop. AiVtot is the upper bound 
microscopic estimate of effective superposition size, i.e., the total number of electrons participating 
in the superposition state. 



Exp. 


Mat. 


v F [ m/ s] 


L[fim] 


si P 




AA tot 


SUNY 


Nb 


1.37- 10 6 


560 


2-3 /iA 


5.5-8.3- 10 9 


3800-5750 


Delft 


Al 


2.02 ■ 10 6 


20 


900 nA 


2.4- 10 6 


42 


UCB 


Al 


2.02 ■ 10 6 


183 


292 nA 


4.23 ■ 10 7 


124 



parameters and corresponding values of AA to t for the three recent experiment [61171 [2Tj 
are listed in table [TJ The SUNY experiment [7J was carried out using a single-junction RF- 
SQUID configuration creating a superposition of two nearly degenerate highly excited states, 
while the Delft \§\ and Berkeley |21| experiments were both made with three-junction flux 
qubits that generated a superposition of degenerate ground states. We also list in table Q] 
the corresponding values for the difference in the two macroscopic observables current and 
magnetic moment between the two branches, respectively SI p and S[i = ASI p , where A is 
the area enclosed by the superconducting loop. 

Electron tunneling through the junction can also contribute to the difference in mode 
occupation in the two branches of the superposition. We have estimated this contribu- 
tion assuming momentum conserving tunneling events |17U18j and find that the tunneling 
dynamics of electrons close to the junction do not significantly alter the superposition size 
estimates given here. We have also analyzed the effects of impurities and defects, i.e., a 'dirty 
superconductor' on the robustness of these flux qubit superpositions. For the high quality 
nanofabricated flux qubit circuits, these effects are restricted to elastic scattering of Cooper 
pairs from impurities. Detailed analysis of the Green's functions given elsewhere [18] shows 
that Sriq is actually independent of the impurity concentration, implying that the results 
obtained here apply both to clean and dirty superconducting flux loops. 

Discussion. — The results in table Q] show that while not trivially small, the num- 
ber of electrons participating in the superpositions is considerably smaller than previous 
estimates that simply counted all electrons within a London penetration depth of the sur- 
face miHIHE]. furthermore, there is a marked contrast between the fact that the values 
of AJV are well below what is generally regarded as macroscopic and the fact that the two 
branches nevertheless have macroscopically distinct values of current and magnetic moment. 
This apparent conundrum arises because the electrons are circulating in opposite directions 
at high speeds (vp) around a path enclosing a relatively large surface area. Our result shows 
that the actual number of electrons that would be found to be behaving differently in the 
two branches if one could measure them at the microscopic level is, however, quite modest. 
This discrepancy between a large difference in the value of an observable quantity (current 
Sip or magnetic moment Sfj,) and the small number of particles actively involved in the 
superposition, ANtot, derives from the indistinguishability and the fermionic character of 
electrons. Thus, given the indistinguishability of electrons, AN to t, be., the total change in 
occupation numbers of all electron modes in the system, is indeed the only meaningful indi- 
cator of how many electrons arc affected when passing from one branch of the superposition 
to the other. In addition, because of the fermion statistics, single electron modes (and also 
Cooper pair modes) corresponding to opposite momenta in the laboratory frame are located 
at diametric extremes on the Fermi surfaces for the two branches. There is a large velocity 
difference between these electrons because vp » v$, even though the number of electrons 
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involved is small. 

We now consider whether the superposition size could be increased to truly macroscopic 
values. One option is to simply increase the number of modes available, by scaling up the 
physical dimensions of the system. Experimental efforts in this direction are underway |22j . 
Another option is to build loops from small cross-section superconducting wires together 
with large area Josephson junctions, for which the larger current differences between the 
superposition branches may be possible |18j . 

The microscopic sizes calculated here for flux superpositions are generally larger than 
typical experimentally realized superposition states in molecular and optical systems [23p4] . 
An intriguing aspect of the numerical estimates is that they are neither trivially small, i.e., 
not unity as expected from a quantum circuit analysis in terms of macroscopic variables , 
nor truly macroscopic. This placement between a microscopic and a macroscopic number 
of participating constituents reflects the key role of the fermi statistics of the electrons that 
emerges from our analysis. Our analysis shows explicitly that for fermionic systems it is 
quite possible to have superpositions which appear large when judged by the magnitude 
of physical observables, but which nevertheless involve only a relatively small number of 
microscopic particles in the superposition. Somewhat surprisingly, this critical role of the 
particle statistics in differentiating these two fundamental properties of large scale superpo- 
sition states has not been appreciated before. Thus, while there is no intrinsic size or number 
scale limiting the existence of macroscopic quantum superpositions, the present work shows 
that the quantum statistics of the constituent particles can nevertheless be important for 
evaluating the effective number of particles participating in a superposition whose branches 
are characterized by macroscopically distinct observables. Whether the corresponding gap 
between difference in observable quantities and number of particles involved is smaller for 
superposition states in a bosonic system is an interesting question for further work. More 
generally, these results imply that existing superconducting flux qubit experiments provide 
realization of quantum superpositions on a scale of a few thousand participating microscopic 
constituent particles. While large in comparison to other experimentally realized superpo- 
sitions, they are not truly macroscopic and thus the formation of superposition states in 
which a macroscopic number of microscopic "elemental" constituents behave differently in 
the two branches remains an open challenge. 
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